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Abstract: The finite temperature physics of the pure glue sector in the improved 
holographic QCD model of |ArXiv:0707.1324| and |ArXiv:0707.T349| is addressed. The 



thermodynamics of 5D dilaton gravity duals to confining gauge theories is analyzed. 
We show that they exhibit a first order Hawking-Page type phase transition. In the 
explicit background of [ArXiv:0707. 13491 , we find = 235 MeV. The temperature 



dependence of various thermodynamic quantities such as the pressure, entropy and 
speed of sound is calculated. The results show a good agreement with the corre- 
sponding lattice data. 
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1. Introduction 

Despite several decades' efforts, an important part of the dynamics of QCD remains 
far from analytical control and in several cases numerical techniques have proved 
too difficult to implement. In particular, recent experiments at RHIC seem to probe 
dynamical properties of the Quark Gluon Plasma (QGP) phase which are not within 
the reach of lattice techniques without extra assumptions. 

On the other hand large- iVc techniques have promised early-on an alternative 
approach to the strongly coupled physics of QCD based on an effective string theory 
description of glue. This route took an interesting twist in 1997 with the advent of 
the Maldacena conjecture [|l|, with the unexpected result that the string theory must 
live in more than four dimensions. In particular there is one extra dimension, known 
as the holographic dimension, that plays the role of (renormalization group) energy 
scale of the strongly coupled gauge theory. 

Since there has been a flurry of attempts to devise such correspondences for 
gauge theories with less supersymmetry with the obvious final goal: QCD. Several 
interesting string duals with a QCD-like low lying spectrum and confining IR physics 
were proposed 0. In the simplest DA example fiavor can be added via the addition 
of Dg branes 0] and its finite temperature phase structure has similarities with QCD 
[|]. Although such theories reproduced the qualitative features of IR QCD dynamics, 
they contain Kaluza-Klein modes, not expected in QCD, with KK masses of the same 
order as the dynamical scale of the gauge theory. Above this scale the theories deviate 
from QCD. Despite the hostile environment of non-critical theory, several attempts 
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have been made to understand holographic physics in lower dimensions in order to 
avoid the KK contamination, based on two- derivative gravitational actions, 

A different and more phenomenological approach was in the meantime developed 
and is now known as AdS/QCD. The original idea described in was successfully 
applied to the meson sector in ^ , and its thermodynamics was analyzed in |^ . The 
bulk gravitational background consists of a slice of AdSs, and a constant dilaton. 
There is a UV and an IR cutoff. The confining IR physics is imposed by bound- 
ary conditions at the IR boundary. This approach, although crude, has been partly 
successful in studying meson physics, despite the fact that the dynamics driving chi- 
ral symmetry breaking must be imposed by hand via IR boundary conditions. Its 
shortcomings however include a glueball spectrum that does not fit well the lattice 
data, the fact that magnetic quarks are confined instead of screened, and asymptotic 
Regge trajectories for glueballs and mesons are quadratic instead of linear. A phe- 
nomenological fix of the last problem was suggested by introducing a soft IR wall, [^]. 
Although this fixes the asymptotic spectrum, it does not allow a proper treatment of 
thermodynamics. In particular, neither dilaton nor metric equations of motion are 
solved. Therefore the "on-shell" action is not really on-shell. The entropy computed 
from the BH horizon does not match the entropy calculated using standard ther- 
modynamics from the free energy computed from the action, etc. Phenomenological 



metrics for the deconfined phase were also suggested, [|I0| |TT| capturing some aspects 
of the expected thermodynamics. 

In [0 an improved model for QCD was proposed. It united inputs from both 
gauge theory and string theory while keeping the simplicity of a two derivative action. 
It could describe both the region of asymptotic freedom as well as the strong IR 
dynamics of QCD. It is a 5d theory like AdS/QCD. 

In this letter we present the finite temperature dynamics in the pure gauge sec- 
tor derived from the setup of We find that this setup describes very well the 



basic features of large-A'^c Yang Mills at finite temperature. It exhibits a first or- 
der deconfining phase transition. The equation of state and speed of sound of the 
high temperature phase are remarkably similar to the corresponding lattice results. 
Moreover, using the zero temperature potential and without adding any extra pa- 
rameter, we obtain a value for the the critical temperature in very good agreement 
with the one computed from the lattice. A detailed derivation of the results will 
appear elsewhere. 



2. Improved Holographic QCD at T=0 

The holographic model introduced in [12] is five- dimensional. The basic fields that 
are non-trivial in the vacuum solution, and describe the pure gauge dynamics, are 
the 5d metric g^^, a scalar $ (the dilaton) that controls the 't Hooft coupling Xt of 
QCD, and an axion a, that is dual to to the QCD 6 angle. Moreover, as the kinetic 
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term of the axion is suppressed by it does not play any role neither in the 

geometry, nor in the evolution of the 't Hooft coupling. It has however a non-trivial 
profile in the vacuum, implying an IR running of the effective ^-angle, ||T^. Quarks 
can be added to the pure gauge theory by adding space-filling D4 — brane pairs 
in the background gauge theory solution. The — D4 tachyon condensation then 
induces chiral symmetry breaking, |TB|, |T2| . 



The action for the 5D Einstein-dilaton theory reads. 



^5 = M^N^ J dJ'x^ R - + V{X) + 2 ^ d^x^ (2.1 



where Mp is the Planck mass ^ and we use the conventions of [|14|. The second term 



in the action is the Gibbons- Hawking with K being the extrinsic curvature on the 
boundary. 

The only nontrivial input in the two- derivative action of the graviton and the 
dilaton is the dilaton potential V{X), where A = e*. A is proportional to the 't 
Hooft coupling of the gauge theory, A = kAj. The constant of proportionality k 
cannot be calculated at present from first principles but as we discuss below all of 
the physical observables turn out to be independent of n. The potential is directly 
related to the gauge theory /^-function once a holographic definition of energy is 
chosen. Although the shape of ^(A) is not fixed without knowledge of the exact 
gauge theory /9-function, its UV and IR asymptotics can be determined. 

In the UV, the input comes from perturbative QCD. We demand asymptotic 
freedom with logarithmic running. This implies in particular that the asymptotic 
UV geometry is that of AdS^ with logarithmic corrections. It requires a (weak- 
coupling) expansion of V^(A) of the form V{X) = 12/^^(1 + t>iA + V2\'^ + ■■■). Here i 
is the AdS radius and Vi are dimensionless parameters of the potential directly related 
to the perturbative /5-function coefficients of QCD, |T^. In conformal coordinates, 
close to the AdS:^ boundary at r = 0, the metric and dilaton behave as 

A„ = -—!— + ... 

log rA 

where the ellipsis represent higher order corrections that arise from second and 
higher-order terms in the /3-function. The mass scale A is an initial condition for the 
dilaton equation and corresponds to Aqcd- 

Demanding confinement of the color charges restricts the large-A asymptotics of 
V{\). In [1^ we focused on potentials such that, as A ^ 00, V^(A) ~ A3 (log A) 



^The physical Planck mass that governs the interactions is MpN^ . We will however call Alp the 
Planck mass for simplicity. 

^We will use a "zero" subscript to indicate quantities evaluated at zero temperature. 
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where a is a positive parameter. The IR asymptotics of the solution in the Einstein 
frame are: 

dsl - e--(i)>r^ + dxl) , Ao e^-^^HTQ ^^""^^ (2.3) 



where the constant C is a positive constant related to A in (|2.2|) . Confinement 
requires a > 1. The parameter a characterizes the large excitation asymptotics of 
the glueball spectrum, ~ ^2(a-i)/a_ p^^, linear confinement, we choose a = 2. 

The parameters of the holographic model a priori are: the Planck mass Mp, 
which governs the scale of interactions between the glueballs in the theory, k, that 
relates A and the 't Hooft coupling, the parameters Vi that specify the shape of 
the potential, the scale A that plays the role of Aqcd and the AdS scale i. The 
latter is not a physical parameter but only a choice of scale: only Ai enters into 



the computation of physical observables. A specific choice for V^(A) was made in [|T2 
with the appropriate asymptotic properties, that only depended on a single parameter 
which can be taken as vi, hence fixing all Vi for i > 1. Furthermore, one can show 
that all of the physical observables both at zero T and finite T are left invariant 
under a rescaling of A. More concretely, given a potential ^(A) and a dilaton profile 
that follows from this potential with an integration constant A, there exists another 
profile with a different integration constant A^ which follows from a rescaled potential 
V^(A) = V{riX) and the two solutions yield the same glueball spectra and the same 
thermodynamic observables. This symmetry allows one to scale away the parameter 
K. Finally, f i and A are fixed by matching to the lattice data for the first two 0"'""'" 
glueball masses. Once A is fixed, all other interesting scales, like the fundamental 
string scale is and the effective QCD string tension a are also fixed. 

This determines all the parameters of the theory except the Planck mass Mpi. 
We shall show below that Mp can be indirectly inferred from the large temperature 
behavior. 



3. The deconfinement transition 

At finite temperature there exist two distinct types of solutions to the action ( |2.1j ) 
with AdS asymptotics, (|2.2|): 

i. The thermal graviton gas, obtained by compactifying the Euclidean time in 
the zero temperature solution with r ~ r + 1/T : 

ds'^ = bl{r) {dr^ + dr^ + dxf) , A = Ao(r). (3.1) 

This solution exists for all T > and corresponds to a confined phase, if the 
gauge theory at zero T confines. 
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ii. The black hole (BH) solutions (in Euclidean time) of the form: 

ds" = h\r) (^-^ + f{r)dr^ + dxl^ , A = A(r). (3.2) 

The function /(r) approaches unity close to the boundary at r = 0. There 
exists a singularity in the interior at r = cxd that is now hidden by a regular 
horizon aX r = where / vanishes. Such solutions correspond to a deconfined 
phase. 

As we discuss below, in confining theories the BH solutions exist only above a certain 
minimum temperature, T > T„iin- 

The thermal gas solution has two parameters: T and A. The black hole solution 
should also have a similar set of parameters: the equations of motion are second order 
for A and /, and first order for h [0. Thus, a priori there are 5 integration constants 
to be specified. A combination of two integration constants of h and A determines A. 
(The other combination can be removed by reparametrization invariance in r). The 
condition / ^ 1 on the boundary removes one integration constant and demanding 
regularity at the horizon, r = r/^, in the form / fh{f^h — r), removes another. The 
remaining integration constant can be taken as fh, related to the temperature by 
AirT = fh. From Einstein's equations one can show [1T3| : 



AnT = b-'{rh){ T ) • (3.3) 



In the large limit, the saddle point of the action is dominated by one of the 
two types of solutions. In order to determine the one with minimum free energy, we 
need to compare the actions evaluated on solutions i. and ii. with equal temperature. 

We introduce a cutoff boundary at r/i = e in order to regulate the infinite 
volume. The difference of the two scale factors is given near the boundary as |ll3| : 



6(e)-6o(e)=C(T)e3 + ... (3.4) 

By the standard rules of AdS/CFT we can relate C{T) to the difference of VEVs of 
the gluon condensate: C(T) oc {TtF'^)t — (TrF^)o. 
The free energy difference is given by 



^ = M^N^^Vi (l5C(T)r^ - nTb^rh)) 

TS 

= 15C(T) M^N^V3r^ - — , (3.5) 

where, in the last equality, we used the fact that the entropy is given by the area 
of the horizon. It is clear that the existence of a non-trivial deconfinement phase 
transition is driven by a non-zero value for the thermal gluon condensate C{T). 

For a general potential we can prove the following statements, that only require 
the validity of the laws of black hole thermodynamics: 
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i. There exists a phase transition at finite T, if and only if the zero-T theory 
confines. 

ii. This transition is of the first order for all of the confining geometries, with a 
single exception described in Hi: 

iii. In the limit confining geometry bo{r) exp(— Cr) (as r oo), the phase 
transition is of the second order and happens at T = 3C/Att. 

iv. All of the non-confining geometries at zero T are always in the black hole phase 
at finite T. They exhibit a second order phase transition at T = 0^ . 



We now sketch a heuristic argument, limited to asymptotics of the type (ZT3). A 



general, coordinate independent proof will appear in [|13 . 

The existence of a minimum black hole temperature T^m in confining theories 
follows from the small and large behavior of the geometries. On one hand, the 
black-hole approaches an AdS-Schwarzschild geometry near the boundary, which 
obeys T = l/nrh- On the other hand, as the horizon approaches the deep interior 
i.e. Vh — > oo, the mass of the black-hole vanishes and the black hole solution ap- 
proaches the zero-T geometry in this limit. In passing, we note that this implies 
vanishing of in this limit. Using the large Vh limit in ( p.3| ), we find the following 
asymptotics for T: 

T-^^r^-\ r.^oo; T ^ — , r, ^ 0. (3.6) 
47r nrh 



The large rh behavior in eq. (|3.6| ) is valid under the assumption that the zero- 
T solution, with IR asymptotics ( |2.3| ), can be continuously deformed into a black 
hole with arbitrarily small mass and arbitrarily large value of r/j. This assumption 
indeed holds, as we will show elsewhere [|I^ for a more general class of confining 
backgrounds. 

Eq. ( p.6| ) shows that for a > 1, that there exists a minimum temperature 
Tmin > above which the black-hole solutions exist. Here, for simplicity, we assume 
a single extremum of the function T(r/j). We illustrate the function T{rh) schemat- 
ically in figure |l]. The simple convex shapes in (a) are due to our assumption of a 
single minimum. In general the function T{rh) may exhibit multiple extrema. Our 



demonstration here can be generalized to these cases |]T3| . In the confining geometries 
a > 1, for a given T > Tmin, there exist a big and a small black hole solution, given 
by rh < rmin and Vh > Tmin respectively, see fig.|l]. The big BH has positive specific 
heat hence it is thermodynamically stable, whereas the small BH is unstable. In the 
borderline confining geometry a = 1, there is a single BH solution. 

Existence of a critical temperature T^. > Tmin for a > 1 follows from the physical 
requirement of positive entropy. From the first law of thermodynamics, it follows that 
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dT/drh = —SdT/drfi. Then, as S* > for any physical system, extrema of J-'{rh) 
should coincide with the extrema of T{r-h). Using also the fact that J-'{rh) — > — oo for 
r/i — > and T{rh) — > near Th oo, we arrive at conclusion (ii) described above: 
There is a first order transition for all of the confining geometries. 

An interesting case is the borderline confining geometry, where coincides with 
Tmin and located at r/^ = oo. The entropy vanishes there because the geometry 
shrinks to zero size. The free energy also vanishes because this point coincides with 
Tc. Therefore the latent heat also vanishes and one has a second order transition. 
Although this geometry is not interesting for the gauge theory, it is of some interest 
for GR. We recall ||12|, that it corresponds to an asymptotically AdS geometry that 
becomes a linear dilaton background in the deep interior. We have shown that such a 
geometry exhibits a second order Hawking-Page transition into a black-hole solution. 
By similar arguments, point iv of the proposition above can also be demonstrated 
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without difficulty. 

Finally, the small Vh asymptotics also allows us to fix the value of the Planck 
mass in (|2.1|) . Small Vh corresponds to high T. This geometry corresponds to an 
ideal gas of gluons with a free energy density JF — > [n'^ /45)N^V3T^. On the other 
hand, as the geometry becomes AdS, eq. ( |3.5| ) implies^ that: J-' 7i^{MpiY N^T'^Vs. 
Hence we conclude that, 

M/ = {A57t^)~K (3.7) 



Using the value of i in [12], we obtain Mp 2.32 GeV. 



4. Numerical Results 

In [|12| an explicit form of the scalar potential with the correct asymptotics was 
proposed. The resulting background, that corresponds to the choice a = 2 in ( p.3| ), 
exhibits asymptotic freedom, linear confinement, and a glueball spectrum in very 
good quantitative agreement with the lattice data. Here we present a numerical 
computation of the relevant thermodynamic quantities in the same theory. 

The potential chosen in was fixed such that the UV expansion reproduces the 
Yang-Mills beta-function up to two loops and has the large-A asymptotics ^(A) ~ 
A^/^(log A)^/^. It depends on two parameters: the first is the overall normalization 
(that fixes the AdS length i and the energy units); the second is bo, that is equivalent 
to the coefficient of linear term in the small A expansion, i.e. Vi. These parameters 
were fit to reproduce the lattice results for the two lowest scalar glueball masses. 

Our general analysis shows that this theory has black hole solutions above a 
temperature T^m and exhibits a first order phase transition at some > Tmin 

To analyze the behavior of the theory at finite temperature, we have solved nu- 
merically Einstein's equations for the metric and dilaton. The integration constants 
were fixed as explained earlier. We find a minimum temperature for the existence of 
black hole solutions, T^m = 210 MeV. 

Next, we compute the free energy difference between the black hole and thermal 
gas solutions, as a function of temperature. As shown in eq. ( p.5|) , there are two 
competing contributions, which must be dealt with separately: 

1. The term 7rTb^{rh) can be obtained directly by evaluating the numerical solu- 
tion at the horizon. 

2. The term 15C{T)i~^ must be extracted by fitting the coefficient of the cubic 
term in the black hole scale factor close to the boundary, b{r) ~bQ{r) ~ C{T)r^. 
This is a large source of error in our numerics, since it is a tiny quantity arising 

difference of 0(1) quantities. 

^It can be shown that the first term in ([3.5D is subleading in the high T hmit. 
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The resulting free energy as a function of the temperature is shown in figure ^, which 
clearly shows the existence of a minimum temperature, and a first order phase tran- 
sition at T = Tc, where J^{T^ = 0. For T < T^, the thermal gas dominates, and the 
system is in the confined phase. For T > T^, the (large) black hole dominates, corre- 
sponding to a deconfined phase. The small black hole branch is thermodynamically 
disfavored at all temperatures. 
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Figure 2: Black hole free energy 
The value we obtain for the critical temperature, Tc = 235 ± 15 MeV, is close 



to the value obtained for large-N Yang-Mills |]T6l, which with our normalization of 
the lightest glueball would be 260 ±11 MeV ^. It should be emphasized that, we did 
not have to adjust any new parameter with respect to the zero-temperature theory 
in order to obtain this result. 

From the free energy we can determine all other quantities by thermodynamic 
identities. However, for numerical precision it is preferable to derive the entropy 
directly as the black hole area, rather than as a derivative of the free energy. The 
latter suffers from the uncertainty in the determination of C{T). Also, due to the 
linear dependence of all thermodynamic quantities on V3, it is convenient to use 
densities. The pressure, and the energy and entropy densities of the deconfined 
phase are given by: 

p = -J^/V3, s = A7iMlNlb%{rn), e = p + Ts. (4.1) 

Next, we present some of the thermodynamic quantities that are compared with 
the lattice results. It is useful to compare dimensionless quantities, so that the 
•^-dependence drops out. 



^The physical units are obtained by fixing mo++ =1475 MeV as in The value 260 ± 11 MeV 
is obtained combining the results in n& and p7| 
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Figure 3: (a) Dimensionless thermodynamic functions and (b) interaction measure. The 
dashed curves correspond to the lattice data of [fl] 



Latent Heat The latent heat per unit volume is defined as the jump in the 
energy at the phase transition, Lh = TcAs(Tc), and it is expected to scale as A''^ in 
the large N^. limit |16|. From eq. ([4.1|) we note that this expectation is reproduced 

0.65^/iVc. This is to be compared 



in our theory. Quantitatively, we find L^^^/T^ 



with the value 0.77 reported in ||16|. 



Equation of state and the interaction measure. A useful indication about 
the thermodynamics of a system is given by the relations between the quantities 
e/T^, 3{p/T^), 3/A{s/T^) (the normalizations are chosen so that they all equal the 
same constant in the case of a free relativistic gas). In figure |^ (a) we compare our 
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results for these quantities with the corresponding lattice results, reported in ||T 



In the low temperature phase, the thermodynamic functions vanish to the leading 
order in N"^ and the jump in e and s at reflects the first order phase transition. 

The interaction measure, (e — 3p)/T'^ (proportional to the trace anomaly), is 
plotted in figure |^ (b), together with the lattice result from [0. From eq. ( |3.5| ), 
e — 3p oc C(T), consistent with our interpretation of C{T) as the gluon condensate. 

Speed of sound. This quantity is defined as = {dp/de)s = s/cy. It is 
expected to be small at the phase transition, and to reach the conformal value = 
1/3 at high temperatures. In figure ^ we compare our results with the lattice data, 
finding good agreement. 
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Figure 4: Comparison between the speed of sound in our model and the lattice result of 
(dashed curves) 



Shear viscosity. In agreement with the general results of |T9[ , the ratio between 
shear viscosity and entropy density is rj/s = (47r)~^. 



5. Discussion 

The model presented here describes well the basic features of laige-Nc Yang Mills 
at finite temperature: it exhibits a first order deconfining phase transition, and 
the temperature dependence of the pressure, entropy, energy density, interaction 
measure and speed of sound in the high temperature phase behave similarly to the 
corresponding lattice results. Without adding any extra parameter, one obtains a 
value for the critical temperature 10 % off the lattice value. 

^These results are for A^^ = 3; we are unaware of similar plots obtained in the large limit. 
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On the other hand the model can be improved in many ways. The latent heat 
Lh/T^ is 40% off the lattice value. Also, our comparison shows that (see e.g. fig 3a) 
approach to the free field limit at high T is slower than the lattice data. This may be 
traced back to the relative smallness of the latent heat in our potential. Although the 
UV and the IR asymptotics of the dilaton potential are fixed by general requirements 
from the field theory, the intermediate region is free to modify. The reason is that the 
low-level glueball spectrum and the thermodynamics near the phase transition are 
not controlled by the same regions of the potential. With a suitable deformation one 
hopes to obtain better agreement with the lattice data. In particular, it is possible 
to obtain a fit to quantities in figs. 3 and 4, well within the errors of the lattice data 



in a temperature range Tc < T < 5Tc 13|. Retrofitting the potential is an interesting 



challenge that we plan to address in |[T3|. 
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